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Recent investigations have looked at the many-body spectra of random two-body interactions. 
In fermion systems, such as the interacting shell model, one finds pairing-like spectra, while in 
boson systems, such as IBM-1, one finds rotational and vibrational spectra. We discuss the search 
for random ensembles of fermion interactions that yield rotational and vibrational spectra, and in 
particular present results from a new ensemble, the "random quadrupole-quadrupole ensemble" 



I. INTRODUCTION: THE PUZZLE OF COLLECTIVITY 

\d '. 

Even-even nuclei have long been known to exhibit a wide range of collective behavior. The broadest classification 
of collective behavior is into three groups: pairing, vibrational, and rotational collectivity[l|i|. The first is in analogy 
with superconductivity in metals, and the latter two based upon quadrupolc deformations of a liquid drop. 
£H ■ The spectral signatures of collectivity are: 

(1) The quantum numbers of the ground state. In particular for pairing one expects JJ S = + . 

(2) Regularities in the excitation energies, in particular among the lowest J T = 2 + ,4 + ,6 + , . . . states relative to the 
t— I ' J n = + ground state. Particularly useful measures are the ratios R42 = E x (4^)/ E x (2^) and i?62 = E x (6i)/ E x (2^). 
C*") , For pairing, vibrational, and rotational collectivity one expects R42 = 1,2, and 3.33, respectively. 

(3) Strong, correlated B(E2) transition strengths, which measure the collectivity of the wavefunctions. (In particular 
one expects strong intraband transitions and weak interband transitions, but we will not consider that further here.) 

■ Traditionally one either compares B(E2) strengths to the single-particle limit (Weisskopf units) or ratios of B(E2)s. 

\ Simple and successful models of collectivity, usually based upon group theory, are well known [2( ; in fact more 

■ sophisticated classifications are possible and are used, but we do not consider them here. 

t-H ' In order to understand the roots of collectivity, start from the coordinate-space Hamiltonian 

"-^ 1 l Kj 

with a chosen nucleon-nucleon potential. Although there are many ways to solve this Hamiltonian, we consider the 
shell- model route: one chooses a single-particle basis {4>i{r)} and writes the Hamiltonian in occupation space using 
second quantization, where a] creates a particle in the single-particle state <pi, 

H = ^ €. L a\a + i ^ Vijkiala^aiak- (2) 

i ijkl 



o 
=5 



- 1—^ 

X 



The ej's are the single-particle energies and the Vy^'s are the two-body matrix elements, in principle computed from 
antisymmeterized integrals of the nucleon-nucleon potential: 

_ 1 

V ijkl = [ d?rcPr'4>*(m{?)V{r, ?) [M^M?) - MVM?)] (3) 



Any interacting shell-model code, such as the one we use, REDSTICK0, reads in: the list of valence single-particle 
states (such as the lsi/2-0<i3/2-0d5/2 or sd space or the IP1/2-I.P3/2-O/5/2-O/7/2 or pf space); information on the many- 
body space, such as valence protons and neutrons; and finally a list of the single-particle energies and two-body matrix 
elements, both of which are read in as numbers. The two-body matrix elements also conserve angular momentum J 
and frequently, but not always, isospin T. 

As mentioned, there exist simple algebraic models of collectivity. Importantly, they have representations both in 
fermion and boson spaces, with the two-body matrix elements derived from the group generators. Furthermore, these 
algebra-based models describe well much of the experimental spectra of even-even nuclides. Therefore, the default 
assumption has long been that the data in turn imply that the real nuclear interaction must have buried deep inside 
strongly "algebraic" components. This assumption can be tested by sampling a large number of Hamiltonians and 
seeing if any signatures of collectivity arise 0. 
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II. THE TWO-BODY RANDOM ENSEMBLE (TBRE) 



The simplest test is to replace the two-body matrix elements with random numbers, independent except that angular 
momentum J (and isospin T) is (are) conserved. This is known as the two-body random ensemble (TBRE) and, as long 
as the ensemble is symmetric about zero, the results are broadly insensitive to the details of the distribution e.g. 
weighting with J, uniform or Gaussian distribution, etc. (The displaced TBRE, or DTBRE 6], which is not symmetric 
about zero, will be discussed briefly below.) The TBRE was originally used to investigate quantum chaosp}, and so 
such investigations considered only states with the same quantum numbers. It was not until a few years ago that 
investigations compared spectral properties of states with different quantum numbers 4] . 

Use of the TBRE is straightforward. A sample interaction from the ensemble is generated, and fed into a shell- 
model diagonalization code. For the TBRE this means that Vjx(ij,kl) — an independent random number. The 
many-body system is referenced by its nuclear physics analog, e.g., " 22 0" mean six "neutrons" (identical fermions) 
in an sd valence space, although the actual many-body space is abstracted and is no longer tied to, for example, the 
radial wavefunctions. The output (for example, the angular momentum of the ground state) is stored in a file for 
later analysis. Typically one compiles several hundred or even several thousand runs to generate good statistics. 

Previous investigations into the properties of the TBRE in fermion systems yielded the following results: 

• Dominance of ground states with J T = + . (Hereafter we will drop parity, as most of the model spaces investigated 
do not have abnormal parity states.) Specifically, in many-body spaces where the fraction of J = states is small, 
typically 4-11 %, after diagonalization of Hamiltonians drawn from the TBRE, the fraction of ground states with 
J = is 45-70%. This is the best known result from random interactions. 

• Pairing-like behavior In addition to J — ground states, one also finds odd-even staggering of ground state 
binding energies, and a ground state "gap." 

For a pairing system, one expects R42 ~ 1. Fig. 1 shows the distribution of R42 over a thousand samples from the 
TBRE for several different many- fermion systems. All are peaked near 1 (pairing), but with broad distributions that 
include R42 = 2 (vibration) and 3.33 (rotational). 

Many-boson systems, in particular the IBM-1, have also been studied using the TBRE, also yielding a dominance 
of J = ground states[jj. In contrast to fermion systems, however, the distribution of R42 is sharply peaked at 2 and 
3.33. Bijker and Frank studied the ratio B(E2 : if -» 2^ )/B(E2 : 2^ -> Of), which turns out for IBM-1 to be strongly 
correlated with R42 and agrees with analytic algebraic limits; plots of such correlations we will call Bijker- Frank plots. 
This is interpreted as evidence for band structure in boso ns sy stems with random interactions. The boson results 
have been successfully explained using a mean-field analysis [lfjj. 

There have been numerous papers "explaining" the fermion results; for reviews see |lll fl2| . Much of the focus has 
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FIG. 2: Correlation of R 62 vs. R42 (a) "realistic" calculations in sd and pf shells; TBRE in (b) 24 Mg; (c) 44 Ti; (d) 4S Ca. 



been on thepredominance of J = ground states; few have considered band structure beyond J = 4, save for a few 
exceptions 0,[l3. In Fig. 2 we give a correlation plot for for i? 62 vs. .R42. Fig. 2(a) shows results using "realistic" 
interactions calculations in the sd 0] and pf shell; 2(b)-(d) are results for the TBRE in several different fermion 
systems. Somewhat surprisingly, we get a strong correlation that follows exactly the naive predictions for pairing, 
vibrational, and rotational collectivity. (Although we do not show it, this result does not appear in single- j systems 
frequently used as testbeds for random interactions). This higher-order band structure is a challenge to any claim to 
understanding random interactions results for fermion systems. Although we do not show it, similar (and sharper) 
results occur in the IBM with random interactions. 

Fig. 3 shows correlations of the ratio B(E2:4 — » 2)/B(E2:2 — > 0) vs. R42 (Bijker-Frank plots), for "realistic" sd and 
pf shell calculations in Fig. 3(a), and for the TBRE in 24 Mg in Fig. 3(b). (For the B(E2) transitions in TBRE, for 
lack of other constraints we assume harmonic-oscillator radial wavefunctions.) While IBM-1 Bijker-Frank plots show 
sharp correlations, the fermion TBRE does not show very strong correlations. (Another approach to collectivity in 
fermion systems with random interactions is the Alaga ratio [lj). 

Fig. 4 shows correlations between B(E2:6 -> 4)/B(E2:2 -> 0) and B(E2:4 -> 2)/B(E2:2 -» 0), which has not been 
previously studied. Again, Fig. 4(a) shows realistic results, while Fig. 4(b) is for the TBRE. The correlations are 
sharper than for the Bijker-Frank plot. Not shown are results for IBM-1, which has still sharper correlations. 

To summarize: fermion systems with TBRE interactions show some signatures of collectivity, but much more weakly 
than boson systems. Why is not understood. This motivates the next section. 



III. THE RANDOM QUADRUPOLE-QUADRUPOLE ENSEMBLE (RQQE) 

One way to probe the difference between random interactions in boson and fermion systems is to search for a 
random ensemble that yield stronger collective structure in fermion systems. One proposal is the displaced TBRE, 
or DTBRE, where the two-body matrix elements are given a constant displacement. While an appealing suggestion, 
we have found that the DTBRE displays collective behavior only for a handful of even-even systems. Even worse, as 
discussed below, the resulting ensemble is not very "random." 

Therefore, inspired by the "consistent-Q" formulation in the interacting boson modelQ, we propose the following. 
Consider the general one-body operator of angular momentum rank 2: 



*i [ & i 



(4) 
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FIG. 3: Bijker-Frank plots for fermion systems, (a) "Realistic" systems in the sd and pf shell, (b)-(d) for 24 Mg: (b) TBRE. 
(c) RQQE with standard quadrupole transition operator, (d) RQQE with 'consistent' transition operator. 

If qij — (i||r 2 Y2||j), then Q is the standard quadrupole operator. But instead we choose the q^ randomly. We then 
define the random quadrupole- quadrupole ensemble to be interactions 

H = -XQ ■ Q. (5) 

The antisymmeterized matrix elements are Vijik = qikqji—quqjk (here and above we have left out details of the angular 
momentum coupling which are straightforward, if tedious, to include). There is also a one-body single-particle energy 
induced by normal ordering of the operators. 

Any member interaction of the RQQE has fewer random parameters than the TBRE. For example, in the sd shell, 
the TBRE has 63 independent, randomly generated two-body matrix elements, while the RQQE, if one assumes 
time-reversal symmetry (that is, q a t = {— l)"' a ~ : "'<fta), has only 5 independent parameters; for the pf shell, the TBRE 
has 195 independent matrix elements and the RQQE has 9. 

The RQQE has J = ground states > 99% for all even-even nuclides we tried. In Figs. 1-4 we show other results 
for the RQQE alongside those for the TBRE. To summarize them: 

• In Fig. 1 we see R42 for the RQQE peaks between 2 and 3.3, that is, somewhere between vibrational and rotational. 

• Although not shown, the correlation of i? 6 2 vs. i? 42 for the RQQE is similar to that of the TBRE in Fig. 2. 

• In Fig. 3(c)-(d), the Bijker-Frank plot shows the correlation between ratios of B(E2)s and i?42- Unlike in the TBRE, 
the RQQE shows much stronger correlation. Fig. 3(c) uses the "standard" quadrupole transition operator, assuming 
harmonic oscillator single-particle states; Fig. 3(d) uses a "consistent" quadrupole transition operator, the same Q 
used in the Hamiltonian. Both results are similar, although the consistent-quadrupole results are sharper. 

• In Fig. 4(c)-(d) we show correlations in B(E2)s in the yrast band 6-4-2-0. The correlations are much sharper than 
for the TBRE, and sharpest for the consistent-quadrupole calculation. 

IV. TEST FOR RANDOMNESS 

Although the RQQE shows stronger spectral signatures for band structure in fermion systems than the TBRE, 
because the former has many fewer random parameters than the latter, we are concerned if the RQQE results are 
truly 'random.' We propose the following as a necessary, if not sufficient, test of randomness. 

For any given many-fermion system, let N be the dimension of the J = subspace. Therefore, any J = 
wavefunction can be represented as a vector of unit length in an iV -dimensional space. If we choose the vectors 
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B(E2:4-> 2)/B(E2:2-> 0) 



FIG. 4: Correlation of B(E2) ratios B(E2: 6 — » 4) / B (E2: 2-^0) versus B(E2: 6 -> 4)/B(E2: 2^0) for (a) "Realistic" 
systems in the sd and pf shell, (b)-(d) for 24 Mg: (b) TBRE. (c) RQQE with standard quadrupole transition operator, (d) 
RQQE with 'consistent' transition operator. 

randomly, then they will uniformly cover a unit sphere, and it is straightforward to show that the angle 9 between 
any two vectors will have a probability oc sin w °~ 2 9. 

We compared ground state vectors for J = ground state. (We ignored interactions that did not have have J = 
ground states. ) We took the dot product between two randomly chosen ground state vectors, computed the angle 
between them, and binned the results, shown in Fig. 5 for 48 Ca. The expected shown for comparison. 

Fig. 5(a) shows the distribution for the fermion TBRE. It follows exactly the expected sin^ 0-2 -' 9 form. 

Fig. 5(b) shows the distribution for RQQE. It is closer to the expected sin^ 0-2 ^ 9 than the DTBRE in (c) and (d) 
but is not completely random. The distortion may be due to the low number of independent random parameters. 

Fig. 5 (c) and (d) shows the distribution for the displaced TBRE (DTBRE) for displacements c ~ 1, 3, respectively. 
(Following the original paper, the width of the displaced Gaussian is 0.6.) Note as displacement c gets larger, the 
distribution moves further from sin^ 0-2 9 and becomes more peaked towards 9 = 0. These ground state wavefunctions 
are not randomly distributed; rather they are clustered about the limit wavefunction for c — > oo. 

This is strong evidence that TBRE wavefunctions are randomly and uniformly distributed; it is also evidence that 
the DTBRE wavefunctions arc not uniformly distributed. The results from RQQE are inconclusive; they are not 
uniformly distributed as the TBRE wavefunctions are, but they are "more random" than DTBRE. 

V. CONCLUSIONS 

We have revisited the question of signatures of collectivity in the spectra of random, two-body interactions in 
fermion systems. Some surprising, previously unknown correlations show up even in the two-body random ensemble 
(TBRE). We proposed a new "random" ensemble, the random quadrupole-quadrupole ensemble (RQQE), and found it 
had significantly sharper signatures of collectivity than the TBRE. While this does not solve the puzzle of collectivity, 
it does demonstrate that typical spectral signatures of collectivity do not rigorously require the standard Elliot 
quadrupole-quadrupole or SU(3) interaction. 

Of course, collective behavior is not surprising if the "ensemble" is not very random, and we looked at the dot 
product between ground state wavefunctions: the RQQE deviates from the expected random distribution but much 
less so than the displaced TBRE or DTBRE. 
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FIG. 5: Distribution of angles between different J = ground state wavefunctions (for 48 Ca) for (a) TBRE; (b) RQQE; (c) 
DTBRE with displacement c = 1; (d) DTBRE with displacement c = 3. Dashed line is sin'^ 0-2 - 1 6, where iVo is the dimension 
of the J = subspace. 
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